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the input required. Sample problems and a complete listing of 
the program is also provided. 
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NOLIN 
A NONLINEAR LAMINATE ANALYSIS PROGRAM 


by 
John J. Kibler 


SUMMARY 


A nonlinear, plane-stress, laminate -analysis program, 
NOLIN, has been developed which accounts for laminae non- 
linearity under inplane shear and transverse extensional 
stress. The program determines the nonlinear stress-strain 
behavior of symmetric laminates subjected to any combination 
of inplane shear and biaxial extensional loadings. The pro- 
gram has the ability to treat different stress-strain be- 
havior in tension and compression, and predicts laminate 
failure uSing any or all of maximum stress, maximum strain, 
and quadratic interaction failure criteria. 

A second program, UNI, has been developed which com- 
putes elastic constants and thermal coefficients of expan- 
Sion for laminae, from constituent properties, to aid in 
compiling input for the NOLIN program. Laminae properties 
can be computed for isotropic or transversely isotropic 
fibers in an isotropic matrix. In addition, the nonlinear 
inplane shear stress-strain curves are computed for the 
laminae by computing the Ramberg-Osgood shear stress para- 
Meter. 

This document provides brief descriptions of both 
programs, a description of the flow of information through 
the NOLIN program, and detailed descriptions of the input 
required for each program. Sections are provided with 
sample problems and sample program output, along with 


complete listings of each program. 


1. NOLIN Program Description 
il: “neroductiron 


The NOLIN program codifies a nonlinear, plane-stress, 
laminate analysis wherein the nonlinear behavior of the 
laminae under inplane shear and transverse extensional 
stresses are taken into account. Both the underlying analy- 
tical development and the computer program are sufficiently 
general to enable the user to study the nonlinear behavior 
of a symmetric laminate subjected to any combination of in- 
plane shear and biaxial extensional loadings. Contained in 
this document are a detailed description of the input re- 
guired to use the NOLIN program, as well as a description of 
the theoretical developments upon which the program is based. 
For a complete theoretical description the reader is referred 
LO Rete ds 

In unidirectional, fiber-reinforced laminae, the trans- 
verse extensional and, particularly, the inplane shear stress- 
strain relationships cannot be accurately characterized as 
linear. The NOLIN program allows for nonlinear representations 
by permitting these stress-strain relationships to take the 
form of Ramberg-Osgood nonlinear relationships. The introduc- 
tion of these nonlinear, Ramberg-Osgood type constitutive 
relationships into a laminate analysis then leads to a set of 
nonlinear equations involving the laminae stress components 
as unknowns. The program then solves this set of equations 
by means of a generalized, Newton-Raphson procedure to give 
the laminae stresses and strains corresponding to the applied 
boundary stresses. 

The theoretical development for this nonlinear, laminate 
analysis incorporates total deformation theory with Ramberg- 
Osgood type stress-strain characterizations to formulate the 
governing nonlinear equations. At the outset the compliance 
tensor is assumed to be the sum of two tensors, the components 


of one are the usual components associated with linear, 
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orthotropic, plane-stress elasticity theory, while the 
second tensor contains the nonlinear elements. By assuming 
a quadratic interaction of the stress components, and re- 
quiring the constitutive relationship to reduce to the re- 
lationships for the uniaxial stress cases of inplane shear 
and transverse extension, the elements of the nonlinear com- 
pliance tensor are explicitly determined. 

Having the nonlinear laminae constitutive relations, 
the usual methods of laminate theory are then utilized to 
obtain the governing, nonlinear equations for the laminate. 

As in linear, laminate theory, the strains of the individual 
laminae are first rotated to a common set of laminate axis, 

and the laminate compatibility relations requiring the corres- 
ponding strains of the individual laminae to be equal are 

then employed. In addition, equilibrium at the laminate 
boundaries is envoked. In this way the required number of 
equations involving the unknown laminae stresses are formulated. 

The program solution procedure for the set of nonlinear 
equations involving the laminae stress components is a Newton- 
Raphson technique generalized to accomodate systems of equa- 
tions. The starting point for the solution procedure is taken 
as the solution of the associated, linear laminate problem, 
where the associated linear problem is obtained by ignoring 
all nonlinear terms. 

Incorporated into the program are three different failure 
criteria, maximum stress, maximum strain and a quadratic inter- 
action criteria. Any or all of these may be employed by the 
program user. Unfortunately, the nonlinear aspects of this 
program preclude the generation of strength envelopes Since 
linear extrapolation is not valid here. Instead a sequence 


of combined loadings may be run. 


The governing equations are formulated so that the 
three stress components in each lamina are the unknowns. Thus, 
for an N-layered laminate, the problem is formulated in terms 
of 3N unknowns. To obtain solutions, 3N equations are then 
required, and these equations consist of three equilibrium 
eguations and 3(N-1) compatibility equations satisfying strain 
compatibility between adjacent laminae. The three equations 


of equilibrium for a laminate under a combined state of stress 


are, 
N 
(k) _ 
TL a eae | 
ot 
N 
(k) = 
i Bigg “Sip = N55 (6) 
ae 
N 
(k) a 
BOyg by FNio 
i: 24 
Where N 


LL’ Noo and Nj 5 are the applied stress resultants, ty 
the thickness of the kth lamina, and subscripts 1 and 2 de- 
note the laminate axes. The 3(N-1) equations of strain 
compatibility are, 
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Equations (6) and (7) are the 3N equations required for the 
solution of the nonlinear laminate problem. When the stress- 
strain relations given by equations (2), (4) and (5) are 
transformed to the laminate reference axes and substituted 


into equations (7), the governing equations can be expressed 


in functional form as, 


F (O54 Snr seer OF, 7 eee) SQ (8) 


1.3 Method of Solution 


Solutions of equations (8) for the 3N stress components 
are obtained by employing a Newton-Raphson iterative scheme. 
The functions B. are first expanded in Taylor series about 
an approximate set of initial stresses, 4. Considering 


only the first order terms of these series, 


By writing, 


where o. are the solution values, equations (9) can be 


rewritten to give 
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For clarity, the notation in equations (10) is, in expanded 
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The solution for se in equation (10) may be taken as the 
approximate, initial stress values for the next iteration 
step, and this process repeated until a result is obtained 
within some desired accuracy. After the stresses are ob- 
tained and transformed to the laminae natural axes, the 
corresponding laminae strains are determined from equations 


(lL), (2) and. (3). 


1.4 Computer Program 


The flow chart. for the computer program is shown in 
Fig. 1. The major sections of the program are the forma- 
tion of the governing equations, the Newton-Raphson solu- 
tion procedure and the failure checks. 

Using the computer program notation, the governing 


equations take the form, 


f[Al- SG + B = SGO (15) 


Where SG and SGO are the stress solution vector and the 
applied stress vector, respectively. A is a matrix of con- 
stant elements which are the coefficients of the linear 
terms in the solution, and B is a vector containing the 
nonlinear terms in the solution. The set of equations (15) 
are equivalent to equations (8). If in equation (15) the 


vector B is set to zero, the resulting equation 


—— 


[A]- SG = SGO (16) 


is the linear laminate solution. The stress vector, SG, 
as determined from equation (16) is taken as the initial 
approximation for the stress vector in the Newton-Raphson 
procedure. 

For the Newton-Raphson procedure it is necessary to 
formulate the derivative of ([A]° SG + B) in equation (15) 


with respect to o. as well as the vector 


DC = (f{A]°*SGO + B - §GO) (17) 


The vector DC corresponds to the vector Fe in equation 
(10), and an explicit evaluation of DC is obtained by 
using the current, approximate value for the solution 
stress vector, SC. The derivative of ([A]-SG + B) is 
designated DB in the computer program, and is equivalent 
to the matrix (AF /30.) os in equation (10). An explicit 
evaluation of DB is also obtained by using the current, 


approximate value for the solution stress vector, SG. 


In the program, the external loading is applied in incre- 
ments. The approximate solution stress vector for the first 
load increment and the first Newton-Raphson iteration is deter- 
mined from equation (16). For the second and third load in- 
crements, the approximate solution stress vectors for the first 
iteration are taken as the final solution stress vectors from 
the previous increments. Solutions for subsequent load in- 


crements are initiated by the following algorithm: 


= * 
(SG; + y)ayrrran = (864) pryar, * (FACTOR) 
(18) 
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The convergence and divergence criteria employed in the 


program are contained in the following expressions: 


MSGiga = SG, )/SG, | ae 


. (19) 
(S65 44 - SG.) /SG, | a | 


Where SG. and SG. are the solution vectors obtained from 
1 1+] th 


the 1 and i+] iterations. Usually values of 
10 ~ and io." are taken for « and A, respectively. However, 
the other values may be input as data to the program. In 
addition, the maximum number of iterations to be allowed is 
input as data. Ten iterations have been found to be suffi- 
cient for most problems. 

The program contains three failure criteria, maximum 
strain, maximum stress, and a quadratic interaction criteria. 
After a solution is obtained for each load increment any or 
all of these failure criterias may be applied to check for 
laminae failure. 

The maximum stress and maximum strain failure criteria 


check, respectively, the laminae stress or strain values in 
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the fiber, and transverse fiber directions against the 
material allowables. These allowables are input to the 


program as data. The quadratic criteria 1s given by 


2 2 2 
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(20) 
tAy > py, Spm + By Spy + Bo Opp .= + 
where the coefficients are functions of the allowable 
stress 
= L 2 -xtllt i 
i P T, C Sa. = - ae 
le” L L 
_ iL ie a 1 
Ago 7 a Boo = 7 2 C (21) 
ened I Aly ol 
al 
A = 
44 (RS)? 
Fr and a are the allowable tension and compression 
stresses in the longitudinal direction, Pe and Fr are the 


allowable tension and compression stresses in the transverse 
direction, and F° is the allowable shear stress. The co- 
efficient Ald is input as data to the program, or a default 
null value is used. 

If a failure criteria is satisfied at the end of a load 
increment, the program determines the failure load through 
linear interpolation. If all failure criteria are being 
checked, and not all indicate failure during the same load 
increment, the program continues loading until all criteria 


indicate failure. 
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Figure 1 - Solution flow Chart 
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2. NOLIN PROGRAM USERS GUIDE 


2.1 Program Description 


This section describes the input data requirements 
for NOLIN (version 2 mod 2). Input procedures have been 
streamlined wherever possible by using one NAMELIST state- 
ment, hence eliminating any input under format control. 

The input flow diagram and input description provide 
all information necessary to specify input data sets cap- 


able of exercising all program options. 
2.2 Input Description 


The initial data required is a message of five cards 
of alphanumeric descriptive information describing the 
problem being solved and printed as a title on the output. 
These five cards may be left blank, but must be included 
ahead of the first NAMELIST deck in the data. This descrip- 
tive message is read only once at the beginning of the pro- 
gram exectuion. The multiple case feature of running 
successive computations is accomplished by supplying multiple 
NAMELIST data sets with the changed variables indicated. 

The following is a description of the input variables 
required for execution of the program. Where appropriate, 
default or suggested values are indicated. The following 
data are supplied through NAMELIST "DATA": 


Program Option Parameters 


TOPT: Ramberg-Osgood parameter sentinel. 
ITOPT=1: Input R-O parameter directly. 
IOPT=2: Determine parameters from stress- 
Strain curve-fit routine. 
COPT: Curve-fit sentinel, (exercised if IOPT=2) 
COPT=1: Same stress-strain data for each 
layer. 
COPT=2: Different sets of data for each 
layer. 13 


EOPT: Exponent option, (exercised if IOPT=2) 
EOPT=1: Determine exponents from curve- 
FILE CrOUtING:. 
EOPT=2: Input exponents to curve-fit 


routine. 


Solution Accuracy Parameters 


KSGM: No. of load incrememts, maximum = 50 
increments. 
SMLT: Load increment multiple, maximum number 


of Newton-Raphson iterations, default 


tS 100: 

EPS: Convergence criteria for Newton-Raphson 
analysis, default value is ares 

UPBD: | Divergence criteria during Newton-Raph- 


son analysis, default value is 20000. 
INMT: Incrementation estimate method, default 


value is 2. 


Layer Description 


NLAY: Number of laminate layers (max. is 20). 
THICK (LAY): Thickness of each layer. 

THETA (LAY) : Orientation of each layer in degrees. 
MATYPE (LAY) : Material kind of each layer (maximum 


number of different materials is 20). 


Material Description 


Ell (MATYPE) : Lamina longitudinal modulus. 

E22 Lamina transverse modulus. 

G12 i Lamina shear modulus. 

VIZ . Lamina major Poisson's ratio. 

S11IT " Lamina longitudinal tensile strength. 
SEL Lamina longitudinal compressive strength. 
S22T - Lamina transverse tensile strength. 


14 


9226 
SiZ 


EPLILT 


EPLIc 
EP22T 


EPZ2C 
GAMA 
Al2 


STY 
sCY 
Jig d 
XM 


XN 


Stress-Strain Data 


LeTS:: 


SLGIIClyMATYPE):: 


pLG22 


SlGlZ 


EPSLL 


EPS22 


EPS12 


(MATYPE) : 


Lamina 
Lamina 
Lamina 


Lamina 


amine 


Lamina 


Lamina 


transverse compressive strength. 
in-plane shear strength. 
longitudinal tensile strain. 
longitudinal compressive strain. 
transverse tensile strain. 
transverse compresSive strain. 


in-plane shear strain. 


Lamina interaction term for quadratic 
interaction criteria - default value 
is 0.0. 


Ramberg-Osgood tension constant 


Ramberg-Osgood compression constant 


Ramberg-Osgood shear constant. 


Ramberg-Osgood shear exponent 
(default value is 3.0). 
Ramberg-Osgood tension exponent 
(default value is 3.0). 


Number 


(input in IOPT=2) 


of stress-strain data points. 


IPTS values of longitudinal lamina 


stresses for each material type. 


TPTS values to transverse lamina 


stresses for each material type. 


IPTS values of in-plane shear stresses 


for each material type. 


IPTS values of longitudinal lamina 


strains for each material type. 


IPTS values of transverse lamina 


strains for each material type. 


IPTS values of in-plane shear strain 


for each material type. 


BS: 


Applied Loading and Failure Criteria 


SOLL: 
SO22: 


SO12: 
TFCN: 


STIFF: 
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Initial axial stress applied to laminate. 
Initial transverse stress applied to 
laminate. 

Initial shear stress applied to laminate. 
Fallure criteria sentinel 

IFCN=1: ultimate stress 

IFCN=2: ultimate strain 

IFCN=3: quadratic interaction 

IFCN=4: all failure criteria 

Ratio of final to initial laminate 
stiffness which constitutes failure 

due to stiffness reduction, default 


value is 0.10. 


3. UNI PROGRAM 


3.1 Description 


Program UNI computes the elastic properties, thermal 
expansion coefficients and the Ramberg-Osgood shear stress 
parameter for the unidirectional fiber bundle or the lamina. 
The fiber may be isotropic or transversely isotropic and the 
matrix is isotropic. The effective elastic properties of 
the composite are calculated from the composite cylinder 
assemblage model proposed by Hashin and Rosen [2] and the 
thermal expansion coefficients from the analytical results 
of Ref. [3]. The program also calculates the Ramberg-Osgood 
shear stress parameter of the matrix as outlined in Ref. [1]. 

The UNI program takes constituent properties as input 
and computes laminae properties as output. The output 
properties provide all the information required as input to 
the NOLIN program. The program accepts families of fiber 
and matrix materials such that an array of laminae proper- 
ties can be generated. This feature can be especially 
useful for sensitivity analyses. All input to UNI is 
accomplished through a single NAMELIST statement "UNID". 


3.2 Input Description 
The following describes the variables required for the 


execution of UNI through NAMELIST "UNID". 


Program Control Variables 


NF: No. of fibers, max. is 20. 
NM: No. of matrices, max. is 20. 
NVM: No. of matrix vol. fractions, max. is 20. 
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Matrix Properties, J=1, NM 


EM(J): 

RHOM (J) : 
ANUM (J) : 
ALPM (J): 


ROMS (J) : 


Young's modulus for Jth matrix. 
Density of Jth matrix. 

Poisson ratio for Jth matrix. 

Coef. of thermal expansion for Jth 
matrix required only if ISTS=l. 

Shear stress Ramberg-Osgood parameter 
for Jth matrix required only if 
NONLIN = 1. 


Isotropic Fiber Properties, J=1, NF 


The following variables are required when ISOT=1: 


Betsy: 

ANUF (J) : 
RHOF (J): 
ALPF (J): 


Young's modulus for dJth fiber. 
Poisson ratio for Jth fiber. 
Density of Jth fiber. 

Coef. of thermal expansion for Jth 


fiber (required only if ISTS=1). 


transversely Isotropic Fiber Properties, J=l NF 
ee, See ee ees Ore Ete ek Ne 


The following variables are required when ISOT=2: 


EFA(J): 
BET (J)s 


ANUPA (J): 


GFA (J): 


ANUFT (J): 


ALPF (J): 


ALPFT (J): 
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Axial Young's modulus for Jth fiber. 
Transverse Young's modulus for Jth fiber. 
Axial Poisson ratio for Jth fiber. 

Axlal shear modulus for Jth fiber. 
Transverse Poisson ratio for Jth fiber. 
Axial thermal expan. coef. for Jth 

Fiber (required only if ISTS=1). 

Transverse: ‘thermal exp... -Coer. "for Jth fiber: 


(required only if ISTS=1). 


Laminae Volume Fractions, J=l1, NVM 


VM(J): Jth volume fraction of matrix material. 


The UNI program computes laminae properties for all 
combinations of fibers, matrix materials, and volume frac- 
tions which are supplied as input. That is, there are a 
total of NF x NM x NVM materials formed from the input. With- 
in a given run the total number of materials (NF xXNMxNVM) 
must be less than 200. 

The following is a description of the properties com- 


puted by UNI and printed as output: 


Calculated Thermo-Elastic Constituent Parameters 


GFT (J): Transverse shear modulus for Jth fiber. 
GF (J): Shear modulus for Jth fiber. 

GM(J): Shear mdoulus for Jth matrix. 

AKF (J) : Plane strain bulk modulus for Jth fiber. 
AKM (J): Plane strain bulk modulus for Jth matrix. 


Fffective Thermo-Elastic Parameters 


AKTS (J): Effective trans. bulk modulus for Jth 
material. 

EAS (J): Effective axial Young's modulus for Jth 
material. 

BUS (id )-* Effective trans. Young's moduius for Jth 
material. 

ANUAS (J): Ff£. Poisson ratio (unidirectional axial 


stress) for Jth material. 

ANUTS (J) : Eff. Poisson ratio (in transverse plane) 
for Jth material. 

GAS (J): Eff. shear modulus (in fiber planes) for 
Jth material. 

GTS (J): Eff. shear mdoulus (in trans. planes) 
for Jth material. 
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ALPAS (J) : 


ALPTS (J) : 


RHOS (J): 
ROCOMP (J): 


Eff. (fiber direction) thermal exp. 

coef. for Jth material. 

Eff. (trans. direction) thermal exp. coef. 
for Jth material. 

Bulk density for Jth material. 
Ramberg-Osgood shear stress parameter 


for Jth material. 


4, SAMPLE PROBLEMS 


The purpose of this section is to present several sample 
problems which illustrate the capabilities of both the UNI 
and the NOLIN programs. Both the program input and output 
are listed at the end of this section to aid in understanding 


the program. 


4.1 UNI - Sample Problems 


eae See Laminae properties for Thornel 50 fibers in two carbon 
matrices have been determined. The transversely isotropic 
fiber option has been used to model the Thornel 50 fibers, and 
a Ramberg-Osgood shear stress parameter is included for the 
two matrix materials. | 
The input data and the resulting output data are given 
in sections 4.1.2 and 4.1.3. Note that the constituent proper- 
ties are mirrored in the program output along with the required 
computed constituent properties. The effective thermoelastic 
properties for the 1-D laminae are printed on the second page 
of output. Units are consistent throughout the program such 
that the units need only be consistent for the constituent 
properties, in this case properties were input as MN/M? . 
The second sample problem combines KEVLAR-49 fibers with 
a range of matrix properties to obtain the 1-D composite prop- 
erties. A list of input data for this case and the correspond-~ 
ing program output are shown in sections 4.1.2 and 4.1.3. It 
is interesting to note that an order off magnitude change in 
matrix modulus results in a factor of five change in trans- 


verse modulus and transverse shear modulus. 


AL. 


4.1.2 UNI SAMPLE PROBLEMS 
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4,2 NOLIN Sample Problems 


4.2.1 The sample problems for the NOLIN program have been chosen 
to exercise several options of the program. Graphite/Epoxy and 
Boron/Aluminum laminates have been modeled, employing material 
properties derived from the Air Force Composites Design Guide. 

The initial problem is a +30° Boron/Aluminum laminate under 
uniaxial tensile loading. Ramberg-Osgood parameters were input, 
with exponents set equal to 3.0, and a maximum strain failure criteria 
was employed. The input data and the program output follow in 
sections 4.2.2 and 4.2.3 respectively. A graph of the axial stress- 
strain curve for the laminate is given in Figure 2. Note that 
the laminate exhibits a non-linear stress-strain behavior from 
loading onset to failure. 

The second sample nonlinear laminate analysis problem is a 
0, +45 Graphite/Epoxy laminate under combined tension and shear 
loading. In this case the Ramberg-Osgood parameters were determined 
by the curve-fit routine from uniaxial stress-strain data supplied 
as input. The input data and program output listing are given 
in sections 4.2.2 and 4.2.3 respectively. Output of the program 
is plotted in the form of axial stress versus axial strain for 


this case in Figure 3. 
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Figure 2 - +30° Boron/Aluminum Laminate - Axial Stress-Strain 
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Figure 3 - Axial Stress-Strain Response of [0+45] Graphite- 
Epoxy Laminate under Axial and Shear Loading. 
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Be Lviw NOLIN SAMPLE PROBLEM INPUT CARDS 


NOLIN DATA 


NOLIN SAMPLE PrmglLems VERS 9 MON 3 


de 0,749,445 LAMINATE, FULL ReOM PAPAM CURPYSeFIT, ALL FAIL CRIT, 
2, *30,°30 LAMTNATF, =O PARAM INPUT, STRATN FAIL CRIT, 
SMATA 
NLAY=3, | 
E11 (1) 2384 ,982sn5, 
E29¢4)2304,90Nk+04, 
VqoCq) eden, 225c4qn, 
Gy oq sse5e7 79k 4nd, 
THICK(4)2942500,59n, 25: 
THETA(1)249,,0,9,-45,, 
lnpts2,FOPT=A1,09pPT=4, 
IPTS=10, 
S$1G61161,1)=6-89F ann, 43, 78F 490, AN, 67F ON, 27, 56F 490,52, 455 490,41, 349F 400, 
48, 236+09,55 (4 26+0%,42 1 E+Nn, FB, INE+IN, 
EPS11C4,1)541, 16-0352, 66-03, 4 60S, 5, 7Ee IS, 7, 6F 9 15,16 209,1, 518-02, 
1,656-)9,2715F-92,2,846-02, 
S1G2204.1) 56 .89F 4,13, 7BF HO, AN BTF HN, 27 SOE $10, Sn FOF HNU 41 S4F HON, 
46, ,236+09,55,47E+09,42 N1E+ONs, F4,9NESDN, | 7 
Erps22(4,1)21,16-93,7 66€203,4,F eS, 5,7E RNG 7 GF aS 1 LE W09,4  SLEMNC, 
1,056-02,2745F 292,2,84EF-02, 
S1G120451)26089F 40,13, 7BEHND AN OTF HN 27 TOFS 10, S41 ADEHNDS 44 S4E FON, 
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NOLIN SAMPLE INPUT CARDS CONT. 


46, 256409, 55 476409, 42,916 +09, FR, GAE+ON, 


oA ie ae en de ee al 2 


5, 3MEMN 28 4 SOF OND,  ARERWD?, 

SAg1e+5 VENI, aNa%ze7, OF HN AST PSA N, 
IfFchs4,1PRINT=1, 

MATYPE(Q) Fish], | 
SaaT1) 84, S2E0S S227 (1087, 9 6Ee01 SI2C1 21, NES 


2yFP447(4)26,68E 90%, 


EPo2T(1)53e77E@N 3, GAMAC1) 31,8276 HN2.S11001 ee, 48 EHO A FSPCOLISA, 79F SNS, 


EP410(1)84,227E eNO FP 20 (4) EL, A5E A027, 
STIFFEQ,1A0e 

Asaolq)s $2, 96273F 2%, 
KSGMs5Q,9MLTF3.9, TTaq 0, INMT SAF 
BNATA 

NLAYs2, 

E4464) 2942,20EF+N95, 
Eeaty,)are1.24bens, 
Vial) e7en ALE SAO, 
Gq2¢1)s2"2,60F+94, 
THICK(4)=2"0,51, 
THETACZIF30 60,7 %H,9, 
IopTsat, 
ST¥(4)s2"1199,Esn9, 
Scy(4)22%446/,c49%,; 
Ty(A)22"93, 08408, 
XMs%,00,XNFS8.0U, 
9741250 1? 

S29 2S 0.5 

SNg?sn,, 

IF CN? 6 


Sq¢1V(t)2l4 00. 0eS22T C1 tO, 7,519 61d HVS T,FRPILT (1220, 75 Ne, 


EPOPT(4)22, 06-027, GAMAC1) 23, 76097, 5911001921430,6 
EP47004q)80,7E-U3,F Pont (4) 27,07, 

STIFFFO.40u, 

KSGMs48,0MLT=4' 
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5. Computer Program Listings 


Source listings of both the UNI and NOLIN computer programs 
follow. The UNI program requires 20 K of computer core storage 
in a CDC 6600 machine, while the NOLIN program requires 60 K of 
core storage. No peripheral devices are required for either 


program for intermediate data storage. 
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